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ABSTRACT 

Last year observations had a profound impact on our views on the amount and 
nature of dark matter in the universe. We give a brief review of the recent 
history of dark matter models beyond the pure cold dark matter universe. In 
view of the most recent cosmological data, we then go on to discuss models with 
a positive cosmological constant. Finally we explicitly analyse a class of particle 
physics models for a dynamical cosmological component with negative pressure 
("quintessence"), in the context of supersymmetric theories. 

1. Introduction 

1.1. The Universe with A 

The year 1998 has witnessed major changes in our perspectives on the amount 
and nature of the different components of dark matter (DM) in the universe. 

The indications for a presently accelerating universe coming from redshift-distance 
measurements of High-Z Supernovae la (SNe Ia)Bi, combined with cosmic microwave 
background (CMB) data! and cluster mass distribution! seem to favour models with 
a non-vanishing cosmological constant A. Indeed the energy density contributed by 
A should roughly be twice as much as the enegy density of the matter of the whole 
universe, thus leading to a flat universe (with fltotai = 1) composed by 1/3 of matter 
and the remaining 2/3 of cosmological constant or vacuum energy. 

If these results have greatly excited the community of astrophysicists and cos- 
mologists, they have not left indifferent the particle physicists. Indeed, most of the 
people belonging to latter group used to think that the famous problem of the small- 
ness of the cosmological constant would have eventually found a solution with the 
discovery of some symmetry that could make it vanishing. On the contrary, the 1998 
observational evidence seems to indicate that A, far from being zero, could instead 
constitute the major source of energy of a flat (critical) universe. 

But what is actually A? Last year data support the view that the universe is 
presently dominated by a smooth component with effective negative pressure; this 
is in fact the most general requirement in order to explain the observed accelerated 
expansion. The most straightforward candidate for that is, of course, a "true" cosmo- 
logical constant!, but a plausible alternative that has recently received a great deal 



of attention is a dynamical vacuum energy given by a scalar field rolling down its 
potential. A cosmological scalar field, depending on its dynamics, can easily fulfill 
the condition of an equation of state w = p/p between —1 (which corresponds to the 
cosmological constant case) and (that is the equation of state of matter). Since it is 
useful to have a short name for the rather long definition of this dynamical vacuum 
energy, we follow the literature in calling it briefly "quintessence" 

The major bulk of this talk will be devoted to the discussion of the "quintessence" 
possibility in the context of supersymmetric (SUSY) theories. But before plunging 
into this, let us remind the readers that the abovementioned data should be taken 
with some caution. 

1.2. The alternatives 

It is indeniably impressive that having three sets of independent data depending 
on two parameters - the fraction densities in matter and cosmological constant, Q m 
and - we get an overconstrained system with a unique solution (what has been 
called with ironical emphasis "cosmic concordance" i). However, one should not un- 
derestimate the fact that for some of the measurements - for instance the position of 
the first acoustic peak in the CMB and the study of the systematics in the SNe la - 
we are still at a preliminary level. While waiting for further observational evidence in 
favour of cosmologies based on cold DM (CDM) and quintessence (also called QCDM 
models), one should then anyway keep an open eye on the alternative solutions to the 
dark matter problem which have emerged after the crisis of the pure CDM standard 
model. 

In the pure CDM model, almost all of the energy density needed to reach the 
critical one (the remaining few percent being given by the baryons) was provided by 
cold dark matter alone. However, some observational facts (in particular the results 
of COBE) put this model into trouble, showing that it cannot correctly reproduce the 
power spectrum of density perturbations at all scales. At the same time it became 
clear that some amount of CDM was needed anyway in order to obtain a successful 
scheme for large scale structure formation. 

A popular option is that of a flat universe realized with the total energy density 
mostly provided by two different matter components, CDM and hot DM (HDM) in 
a convenient fraction. These models, which have been called mixed DM (MDM)l, 
succeeded to fit the entire power spectrum quite well, although - having Q m = VL tot = 1 
- they obviously can't account for the most recent data. 

Another interesting possibility for improving CDM models consists in the intro- 
duction of some late time decaying particle!. The injection of non-thermal radiation 
due to such decays and the consequent increase of the horizon length at the equiv- 
alence time could lead to a convenient suppression of the excessive power at small 
scales (hence curing the major desease of the pure CDM standard model). As ap- 



pealing as this proposal may be from the cosmological point of view, its concrete 
realization in particle physics models meets several difficulties. Indeed, after consid- 
ering cosmological and astrophysical bounds on such late decays, it turns out that 
only few candidates survive as viable solutions (for a recent analysis in the context 
of SUSY extensions of the SM with or without R parity see Refill). 

Another route which has been followed in the attempt to go beyond the pure 
CDM proposal is the possibility of having some form of warm DM (WDM). The 
implementation of this idea is quite attractive in SUSY models where the breaking 
of SUSY is conveyed by gauge interactions instead of gravity (these are the so-called 
gauge mediated SUSY breaking (GMSB) models). In these schemes the gravitino 
mass loses its role of fixing the typical size of soft breaking terms and we expect it to 
be much smaller than in the more traditional supergravity models. The gravitino in 
GMSB theories can behave as a WDM candidate. Unfortunately, critical universes 
with pure WDM are known to suffer from serious troublesS. It has been shown 
that even variants of the pure WDM models (with an additional HDM component or 
with a non-vanishing cosmological constant) still exhibit various problems in correctly 
reproducing the large scale structure datac3. 

Finally, the alternative to pure CDM which succeeded to fit the whole power 
spectrum as well as the MDM models was provided by the so-called xCDM models^. 
In those schemes, CDM is accompanied with a cosmological constant-like contribution 
(named x) adding up to reach the critical energy density. In view of the recent 
observations that we presented above, it is clear that a;CDM models today enjoy 



the ma 



or success. We now turn to discuss these schemes in the context of SUSY 



theories il. 

2. Supersymmetry and Quintessence 

As already said in the introduction, the 1998 indications for an accelerating uni- 
verse have recently drawn a great deal of attention on cosmological models with 
Q m ~ 1/3 and ~ 2/3. More generally, the role of the cosmological constant in 
accelerating the universe expansion could be played by any smooth component with 
negative equation of state Pq/pq = wq ^ — 0.6i@, as in the so-called "quintessence" 
models (QCDM)i, otherwise known clS X CDM mo delsi. 

2.1. Scalar field cosmology 

A natural candidate for quintessence is given by a rolling scalar field Q with 
potential V(Q) and equation of state 

Q 2 /2 - V(Q) 
WQ Q 2 /2 + U(Q)' 

which - depending on the amount of kinetic energy - could in principle take any 



value from —1 to +1. The study of scalar field cosmologies has shownElO that for 
certain potentials there exist attractor solutions that can be of the "scaling" or 
"tracker" Bel type; that means that for a wide range of initial conditions the scalar 
field will rapidly join a well defined late time behavior. 

If pq <C pb, where ps is the energy density of the dominant background (radiation 
or matter), the attractor can be studied analytically. 

In the case of an exponential potential, V(Q) ~ exp (— Q) the solution Q ~ hit is, 
under very general conditions, a "scaling" attractor in phase space characterized by 
Pq/Pb ~ const@'0H This could potentially solve the so called "cosmic coincidence" 
problem, providing a dynamical explanation for the order of magnitude equality be- 
tween matter and scalar field energy today. Unfortunately, the equation of state for 
this attractor is wq = wb, which cannot explain the acceleration of the universe 
neither during RD {w ra d = 1/3) nor during MD {w m = 0). Moreover, Big Bang 
nucleosynthesis constrain the field energy density to values much smaller than the 
required ~ 2/3 0£H. 

If instead an inverse power-law potential is considered, V(Q) = M A+a Q~ a , with 
a > 0, the attractor solution is Q ~ where n = 3(wq + 1), m — S(wb + 1); 

and the equation of state turns out to be wq = (wb«-2)/(« + 2), which is always 
negative during MD. The ratio of the energies is no longer constant but scales as 
Pq/Pb ~ a m ~ n thus growing during the cosmological evolution, since n < m. pq 
could then have been safely small during nucleosynthesis and have grown lately up to 
the phenomenologically interesting values. These solutions are then good candidates 
for quintessence and have been denominated "tracker" in the literature 0'0H. 

The inverse power-law potential does not improve the cosmic coincidence problem 
with respect to the cosmological constant case. Indeed, the scale M has to be fixed 
from the requirement that the scalar energy density today is exactly what is needed. 
This corresponds to choosing the desired tracker path. An important difference exists 
in this case though. The initial conditions for the physical variable pq can vary 
between the present critical energy density p° cr and the background energy density ps 
at the time of beginningil (this range can span many tens of orders of magnitude, 
depending on the initial time), and will anyway end on the tracker path before the 
present epoch, due to the presence of an attractor in phase spaceSEl. On the contrary, 
in the cosmological constant case, the physical variable p\ is fixed once for all at the 
beginning. This allows us to say that in the quintessence case the fine-tuning issue, 
even if still far from solved, is at least weakened. 

A great effort has recently been devoted to find ways to constrain such models 
with present and future cosmological data in order to distinguish quintessence from 
A models^H. An even more ambitious goal is the partial reconstruction of the scalar 
field potential from measuring the variation of the equation of state with increasing 
redshiftil. 

On the other hand, the investigation of quintessence models from the particle 



physics point of view is just in a preliminary stage and a realistic model is still 
missing (see for example Refs. 001111). There are two classes of problems; the 
construction of a field theory model with the required scalar potential and the inter- 
action of the quintessence field with the standard model (SM) fields^. The former 
problem was already considered by Binetruyil, who pointed out that scalar inverse 
power law potentials appear in supersymmetric QCD decries (SQCD)H with N c col- 
ors and Nf < N c flavors. The latter seems the toughest. Indeed the quintessence field 
today has typically a mass of order H ~ 10~ 33 eV. Then, in general, it would medi- 
ate long range interactions of gravitational strength, which are phenomenologically 
unacceptable. 

In the remaining part of the talk, we will address in more details these problems 
in the framework of SQCD, following the work done in Ref.0. 



2.2. Susy QCD 

As already noted by Binetruy@, supersymmetric QCD theories with N c colors and 
Nf < N c flavors!! may give an explicit realization of a model for quintessence with 
an inverse power law scalar potential. The remarkable feature of these theories is 
that the superpotential is exactly known non-perturbatively. Moreover, in the range 
of field values that will be relevant for our purposes (see below) quantum corrections 
to the Kahler potential are under control. As a consequence, we can study the 
scalar potential and the field equations of motion of the full quantum theory, without 
limiting ourselves to the classical approximation. 

The matter content of the theory is given by the chiral superfields Qi and Q i 
(i = 1 . . .Nf) transforming according to the N c and N c representations of SU(N C ), 
respectively. In the following, the same symbols will be used for the superfields Qi, 
Qi, and their scalar components. 

Supersymmetry and anomaly-free global symmetries constrain the superpotential 
to the unique exact form 

/ A (3N C -N f )\ N c -N f 

* = (».-»i)[-toar-) W 

where the gauge-invariant matrix superfield = Qi ■ Qj appears. A is the only 
mass scale of the theory. It is the supersymmetric analogue of A.qcd, the renor- 
malization group invariant scale at which the gauge coupling of SU(N C ) becomes 
non-perturbative. As long as scalar field values Qi,Qi ^> A are considered, the the- 
ory is in the weak coupling regime and the canonical form for the Kahler potential 
may be assumed. The scalar and fermion matter fields have then canonical kinetic 
terms, and the scalar potential is given by 

_ Nf 1 
V(Q h Qi) = £ (\F Q f + \F^f) + -D a D a (2) 



where Fq. = dW/dQi, Fq = OW/dQ^ and 



D a = Q\t a Q t - Q^Ql . (3) 

The relevant dynamics of the field expectation values takes place along directions in 
field space in which the above D-term vanish, i.e. the perturbatively flat directions 
(Qia) = (Qia)i where a = 1 • • ■ N c is the gauge index. At the non-perturbative level 
these directions get a non vanishing potential from the F-terms in which are 
zero at any order in perturbation theory Gauge and flavor rotations can be used to 
diagonalize the (Qi a ) and put them in the form 

m \ - /7T f \ - qiSia 1 - a - N f 
Wia) WiJ o N f <a<N c - 

Along these directions, the scalar potential is given by 

_ A 2a /N f 

= (V(Q u Q i )) = 2 „ t ■ - E 



2 



3N C -N f 1 

a = -, a = . 

N c - N f ' N c — N f 

In the following, we will be interested in the cosmological evolution of the Nf expec- 
tation values qi, given by 

dV 

(Q i + 3HQ t + —r) = , i=l,...,JV>. 
dQl 

In RefE the same initial conditions for all the Nf VEV's and their time derivatives 
were chosen. With this very peculiar choice the evolution of the system may be 
described by a single VEV q (which we take real) with equation of motion 

thus reproducing exactly the case of a single scalar field $ in the potential V = 
A 4+2 9$ -2 9/2 considered in Refsffli. In this paper we will consider the more general 
case in which different initial conditions are assigned to different VEV's, and the 
system is described by Nf coupled differential equations. Taking for illustration the 
case Nf = 2, we will have to solve the equations 



with H 2 = 8ii/3Mp (p m + p r + Pq), where Mp is the Planck mass, p m ( r ) is the matter 
(radiation) energy density, and pq = 2{q 2 + gf) + f(gi, ^2) is the total field energy. 



2.3. The tracker solution 

In analogy with the one-scalar case, we look for power-law solutions of the form 

g^ = CW Pl , i = l,---,N f . (6) 

It is straightforward to verify that - when pq <C p B - the only solution of this type 
is given by 



1 — r 

Pi=P= — , Q = C 

with 

X 



X l-r p2(3-r) 

4 m (1 + r) 



1/4 



N f , 



[1 -r) 2 [12-m(l + r)] ' 

where we have defined r = Nf/N c (= 1/N C , . . . , 1 — 1/N C ). This solution is charac- 
terized by an equation of state 

1 + r 1 - r 
W Q = -^T w b if~ ' ( ' 

Eq. (0) can be derived as usual from energy conservation i.e. d (a 3 pg) = —3 a 2 pQ. 

Following the same methods employed in Ref.0 one can show that the above 
solution is the unique stable attractor in the space of solutions of eqs. (||). Then, even 
if the qi's start with different initial conditions, there is a region in field configuration 
space such that the system evolves towards the equal fields solutions (B) , and the 
late-time behavior is indistinguishable from the case considered in Ref.c3. 

The field energy density grows with respect to the matter energy density as 

PQ MI+l> , Q , 
— ~ a 2 ( (8) 

Pm 

where a is the scale factor of the universe. The scalar field energy will then eventually 
dominate and the approximations leading to the scaling solution will drop, so that 
a numerical treatment of the field equations is mandatory in order to describe the 
phenomenologically relevant late-time behavior. 

The scale A can be fixed requiring that the scalar fields are starting to dominate 
the energy density of the universe today and that both have already reached the 
tracking behavior. The two conditions are realized if 

v(qo) - Plu , v"(q ) ~ H 2 , (9) 



where p° crit = 3MpHq/8it and q are the present critical density and scalar fields VEV 
respectively. Eqs. (Dl) imply 
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(10) 
(11) 



Depending on the values for Nf and iV c , A and (fa /A assume widely different values. 
A takes its lowest possible values in the N c — > oo (JV/ fixed) limit, where it equals 
4 • lQ-^tf/N}) 1 ' 6 GeV (we have used p° crit / M p = ( 2 -5 • lO- 31 /* 1 / 2 ) 4 ). For fixed V c , 
instead, A increases as Nf goes from 1 to its maximum allowed value, Nf = 1 — iV c . 
For N c ^ 20 and Nf close to N c , the scale A exceeds Mp. 

The accuracy of the determination of A given in fllPf) depends on the present error 
on the measurements of H , i.e., typically, <5A/A = |^ 5H /H ^ 0.1. 

In deriving the scalar potential (0) and the field equations @ we have assumed 
that the system is in the weakly coupled regime, so that the canonical form for 
the Kahler potential may be considered as a good approximation. This condition is 
satisfied as long as the fields' VEVs are much larger than the non-perturbative scale 
A. From eqs. (|T0| ) and (|TTD, one can compute the ratio between the VEVs today and 
A, and see that it is greater than unity for any Nf as long as N c ^ 20. 



2.4- Interaction with other cosmological fields 

The superfields Qi and Q t have been taken as singlets under the SM gauge group. 
Therefore, they may interact with the visible sector only gravitationally, i.e. via non- 
renormalizable operators suppressed by inverse powers of the Planck mass, of the 
form 

QiQi 



d A e kh 



Ml 



where 



(12) 



1TT|) we know that 



bj represents a generic standard model superfield. From 
today the VEVs are typically 0(M P ), so there is no reason to limit ourselves to 
the contributions of lowest order in |Q| 2 /Mp. Rather, we have to consider the full 
(unknown) functions (3^ and the analogous /3's for the Q^s. Moreover, the require- 
ment that the scalar fields are on the tracking solution today, Eqs. @ implies that 
their mass is of order ~ Hq ~ 10~ 33 eV. 

The exchange of very light fields gives rise to long-range forces which are con- 
strained by tests on the equivalence principle, whereas the time dependence of the 
VEVs induces a time variation of the SM coupling const ant sSH. These kind of 
considerations sets stringent bounds on the first derivatives of the 
today, 



s and f3 Jh s 



■Jt 



a 



3 1 
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where Xi = qi/Mp. To give an example, the best bound on the time variation of 
the fine structure constant comes from the Oklo natural reactor. It implies that 
\a/a\ < 1CT 15 yr -1 (see Ref.il), leading to the following constraint on the coupling 
with the kinetic terms of the electromagnetic vector superfield V, 



Vi a Vi < 10" 6 



Mr 



(13) 



where (<ji) is the average rate of change of q^ in the past 2 x 10 9 yr. 

Similar -although generally less stringent- bounds can be analogously obtained 
for the coupling with the other standard model superfieldsil. Therefore, in order to 
be phenomenologically viable, any quintessence model should postulate that all the 



unknown couplings /3 Jl 's and 1 's have a common minimum close to the actual value 
of the qi's? 

The simplest way to realize this condition would be via the least coupling principle 
introduced by Damour and Polyakov for the massless superstring dilaton in Ref£3, 
where a universal coupling between the dilaton and the SM fields was postulated. In 
the present context, we will invoke a similar principle, by postulating that ft 1 = f3 
and f3 31 = f3 for any SM field 4>j an d any flavor i. For simplicity, we will further 
assume (3 = (3 . 

The decoupling from the visible sector implied by bounds like (0) does not neces- 
sarily mean that the interactions between the quintessence sector and the visible one 
have always been phenomenologically irrelevant. Indeed, during radiation domination 
the VEVs qi were typically Mp and then very far from the postulated minimum 
of the /3's. This leads, for the quintessence fields, to the generation of SUSY breaking 
masses (proportional to H) by the same mechanism discussed by Dine, Randall, and 
Thomas in Ref.H. 

The main phenomenological effect of these time- dependent SUSY breaking masses 
is to prevent the fields from taking too large values. This results in an improved 
attraction to the tracker solution (for further details on this point see Ref.0). 



2.5. Numerical Results 



The general results of the previous discussion are illustrated in the figures for the 
particular case Nf = 2, N c = 6. 

° An alternative way±o suppress long-range interactions, based on an approximate global symmetry, 
was proposed in Ref£3. 
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Fig. 1. The evolution of the energy densities p of different cosmological components is given as a 
funcion of red-shift. All the energy densities are normalized to the present critical energy density 
p° r . Radiation and matter energy densities are represented by the short-dashed lines, whereas the 
solid line is the energy density of the tracker solution discussed in Section 2.3. The long-dashed line 
is the evolution of the scalar field energy density for a solution that reaches the tracker before the 
present epoch; while the dash-dotted line represents the evolution for a solution that overshoots the 
tracker to such an extent that it has not yet had enough time to re-join the attractor. 

In Fig.l the energy densities vs. redshift are given. We have chosen the same 
initial conditions for the two VEVs, in order to effectively reproduce the one-scalar 
case of eq. ([|), already studied in Refs.0'0i!. No interaction with other fields of the 
type discussed in the previous section has been considered. 

We see that, depending on the initial energy density of the scalar fields, the tracker 
solution may (long dashed line) or may not (dash-dotted line) be reached before the 
present epoch. The latter case corresponds to the overshoot solutions discussed in 
Ref.@, in which the initial scalar field energy is larger than ps and the fields are 
rapidly pushed to very large values. The undershoot region, in which the energy 
density is always lower than the tracker one, corresponds to p Q cr < Pq < p™. Thus, 
all together, there are around 35 orders of magnitude in p l q at redshift z + 1 — 10 10 for 
which the tracker solution is reached before today. Cleary, the more we go backwards 
in time, the larger is the allowed initial conditions range. 

Next, we explore to which extent the two-field system that we are considering 
behaves as a one scalar model with inverse power-law potential. In Fig. 2 we plot 
solutions with the same initial energy density but different ratios between the initial 
values of the two scalar fields. Given any initial energy density such that - for 
QT/lT = 1 ~~ the tracker is joined before today, there always exists a limiting value 
for the fields' difference above which the attractor is not reached in time. 
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Fig. 2. The effect of taking different initial conditions for the fields, at the same initial total field 
energy. Starting with q" L /q™ — 10 15 the tracker behaviour is not reached today. For comparison 
we plot the solution for q^/q™ = 1- 
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